Consistent Thermodynamics for Quasiparticle Boson System with 

Zero Chemical Potential 

Shaoyu Yin^* and Ru-Keng Su^'^'^'l" 

1. Department of Physics, Fudan University, Shanghai 200433, P. R. China 

2. CCASTfWorld Laboratory), P.O.Box 8730, Beijing 100080, P. R. China 

3. Center of Theoretical Nuclear Physics, 
National Laboratory of Heavy Ion Collisions, Lanzhou 730000, P. R. China 

Abstract 

The thermodynamic consistency of quasiparticle boson system with effective mass m* and zero 
chemical potential is studied. We take the quasiparticle gluon plasma model as a toy model. The 
failure of previous treatments based on traditional partial derivative is addressed. We show that 
a consistent thermodynamic treatment can be applied to such boson system provided that a new 
degree of freedom m* is introduced in the partial derivative calculation. A pressure modification 
term different from the vacuum contribution is derived based on the new independent variable m* . 
A complete and self-consistent thermodynamic treatment for quasiparticle system, which can be 
widely applied to effective mass models, has been constructed. 
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I. INTRODUCTION 



System with complicated interaction is very common in real world and rich in physical 
phenomenon but very hard to deal with. In practise, physicists have introduced many theo- 
retical tools, such as the Hartree-Fock approximation, approximate secondary quantization, 
summation of Feymann diagrams, etc., among which the quasiparticle approximation is one 
of the simplest first order approximation. Most of the quasiparticle approximations take the 
system as consisting of quasiparticles without interaction, while the effect of interaction is 
represented by the effective mass m*, usually being a function of temperature T and parti- 
cle density p. It should be emphasized that the effective mass quasiparticle model is only 
a rough approximation since many detailed properties and higher order contributions are 
neglected. However, to some degree, it is still an efficient tool to discuss some macro prop- 
erties of complex system because of its simplicity. The application of quasiparticle model to 
;he study of thermodynamics can be found in many references, such as to the quark system 



ii,y,y,y,y,y,0,y,y,Q 



or quark-gluon plasma (QGP) [U 



12, Q, 14, H, 



16 



17 



18|. 



However, it is well known that the effective mass quasiparticle model suffers from the 
thermodynamic inconsistency. To illustrate this inconsistency, for a quasiparticle boson sys- 
tem with chemical potential /i = 0, the standard expressions of pressure, p = —Q/V = 



I Tj^T-^dk, will 

exp(e/T)— 1 ' 

. The problem comes from 



— ^ /o^ ^"^[^ ~ exp(— e/T)](iA;, and internal energy density, e = /q°° k 
not satisfy some differential relations, such as e = — p 
that the effective mass, being a function of medium parameters, entangles the usual ther- 
modynamic relations via the relativistic dispersion relation e(A;) 



'/c^ + m*2(T, p). As was 
shown in Ref.[l9j, the usual thermodynamic treatment with traditional partial derivative to 
obtain the thermodynamic quantities is problematic. The reason is that the thermodynamic 
potential f2 of quasiparticle model is not only a function of characteristic variables T, V and 
/i, but also depends explicitly on the effective mass, = f2(T, V, /i, m*). To overcome this 
difficulty, there are many different treatments in the market. For example: (I). Calculate 
the derivatives of m* with respect to T and p in the differential relations for a reversible 
process, 

(M = -SdT - pdV - Ndp, (1) 



on' 

df 



p 



dp) 



(2) 



T,V 



2 



where S, p and are entropy, pressure and average particle number, respectively. Many 

an dm* on dm* 

dm* dp dm* dT 



extra terms involving and ^Sr will emerge in the expressions of p, e and entropy 



191 , though many authors follow this 



density s |J, Q, Q. But as was summarized in Ref. 
direction, the extra terms in different references contradict each other. (II). Introduce a 
temperature- and/or density-dependent vacuum energy [g], [is], 12\ and force the vacuum to 



cancel the thermodynamic inconsistency. However, the physics of the vacuum is still unclear. 



(III). A new treatment had been suggested in our previous paper [19|. We introduce a new 
degree of freedom m* for the quasiparticle model effectively in the thermodynamic derivative 
relation and rewrite Eq.(l) as 



dn = -SdT -pdV - Ndfi + Xdm*, (3) 



then Eq.(2) becomes 



(4) 

where X is an extensive quantity corresponding to the intensive variable m* . We emphasize 
that after introducing the term Xdm* in Eq.(3), we have Eq.(4) instead of Eq.(2). Since 



m* is an invariant in the partial derivatives for S, p and A^, a. 



1 the extra terms involving 



and ^^^fr are forbidden there. We have proven in Ref. 1^, by keeping m* invariant 
during the derivatives, the results remain self-consistent and in accord with the equilibrium 
statistics. 

Noticing that at a fixed instant in the reversible process, the system is in an equilibrium 
state with temperature and density denoted as To and po? respectively, then the effective 
mass of the quasiparticle becomes constant m*(To,po) = In this equilibrium state, the 
system reduces to a usual ideal gas system with constant particle mass mo, whose thermody- 



namic quantities can be directly obtained statistically. In Ref.[l9(] we have proven that the 
thermodynamic quantities calculated in equilibrium state are just the corresponding quan- 
tities given by Eqs.(3) and (4). The introduction of the new variable m* can be understood 
in the following aspect. In the quasiparticle approximation, the effective mass summarizes 
the interaction, the confinement mechanism, etc., which must be reflected in a consistent 
thermodynamic treatment. It forces us to introduce m* as a new independent variable 
in the partial derivative to get a suitable treatment. Employing the quark mass density- 



3 



dependent model, we had shown this treatment is thermodynamically self-consistent 19|. 
The ambiguities and inconsistencies in previous thermodynamic treatments are overcome. 

The motivation of this paper is to extend our study from the fermion system to the 
boson system with non-conserved particles. The chemical potential fi equals zero for such 
system. We here take a quasiparticle gluon plasma (qGP) model, which was first introduced 
m Refs.[ll|, Q, toy model to illustrate our treatment. This model is quite rough: 
On the one hand, it can not describe the properties of the real QGP exhibited by recent 



RHIC exp eriments 20|, l2l|]. One the other hand, one can calculate the hard thermal or dense 



loops [22, l23|, I2J] by temperature quantum theory and get higher order contributions beyond 



mean field and quasiparticle picture. However, we will still employ this model. The reasons 
are threefold: (1). The model is very simple and we can use this model to illustrate our 



treatment transparently. (2). This model has been emp 



thermodynamic inconsistency in the physics market 



13 



oyed by many authors to discuss the 



17 



18l | . We use this model in favor 



of c omp aring our treatment with others. (3). This model has the lattice simulation results 



25 



261], which can be used as standards in comparing different treatments to expose their 



advantages and shortcomings. We will focus our attention to treatment (I) and treatment 
(III) since the study of qGP model by treatment (II) was shown in Ref. [isl in detail. We will 
prove that the results given by traditional partial derivative treatment (I) are completely 
different from those of the lattice simulation. This treatment is failure for the qGP system 
as for the quark system jigj ]. 

Another important problem involves the negative pressure. Of course, a real gluon plasma 
has no negative pressure. But as is predicted by lattice simulation, the QGP system ex- 
periences a phase transition at a critical temperature below which the system becomes 
comfined. One need to introduce the negative pressure at T < Tc indicating the instab ility 
in the quasiparticle model to describe the phase transition at Tc. In treatment (I) 0, [?, 10], 



the authors claimed that the term -^^^w- contributes a negative pressure which can cause 

am* op I 1 ° ^ 

the instability at low T. But unluckily, in Ref. 19[, we proved this term can not exist in 
a consistent treatment. In treatment (II), the authors obtained the negative pressure from 
the vacuum , \v\ . In this paper, we will show that the negative pressure can exist nat- 
urally in treatment (HI). It is not necessary to involve the vacuum. The newly independent 
variable m* permits a freedom for an arbitrary function f{m*) in the formula of pressure, 
whose value can be determined by the comparison with experiment results or lattice data. 
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The paper is organized as follows: In the next section, using the lattice data as the 
criterion for comparison, we will examine the qGP model by treatment (I). In Sec. Ill, similar 
comparison will be made between the lattice data and our treatment (III) . Then in Sec. IV, 
a modification term to pressure and entropy is derived within our self-consistent treatment 
(III). We will show this modification term is different from the vacuum contribution but a 
natural result of our newly introduced variable m*. Sec. V is a brief summary. 



II. THE RESULTS OF TREATMENT (I) 

For the thermodynamic system with effective mass, after calculating the derivatives of 
m* to T and p in Eq.(2), one obtained [lO] 

_ ^ yd^ dm* 
^ dV ^ dm* dp 

~ an an ^ on dm* 

dQ dQ dm* 

where VL = VL/V denotes the thermodynamic potential density. The qGP model is a perfect 
boson system with vanishing chemical potential p and temperature- and density-dependent 
particle mass m*, so the thermodynamic potential density 

gT 



Q = y In (l - e-^^^^^) dk\ (8) 

which is explicitly independent of V for infinite qGP system. Eqs.(5) and (6) reduce to 

— _^ '^^* (9) 

P P Qjyi* Qp ' 

, = n_Tg-T#^. (10) 

dT dm* dT ^ ' 

Now we employ Eqs.(7-10) and three different models to study the properties of qGP model. 
The results are shown in Figs. 1-3. 

The first model is an SU(2) qGP model. The detail of this model has been addressed in 

n 

Ref. [13| and all the parameters including the fitting of m* which we used for calculation are 
from there. The results are shown in Fig.l, where the curves of energy density, pressure and 
entropy density are normalized to their Steven-Boltzmann limits, respectively. The solid 



curves adopted from Ref. |25l| are given by the lattice QCD calculation which we use as the 
standard to compare different treatments. The dotted curves are the results of Eqs.(7-10). 
We see from Fig.l that the deviation between solid curves and the corresponding dotted 
curves are remarkable, especially in low temperature regions. 

The second and the third models are SU(3) qGP models. The details of these models 



have been addressed in Ref.flSj, where the effective mass of SU(3) gluon is given by 

*2 Qoas(T)p 

m = — , model I, (11) 

and 



^2 3aoas(T)p 
m = 



2^ , modem, (12) 

where the constant = 2.15 and the two-loop order running coupling constant 

27r 511n(21n(T/A^)) ^^ 
' llln(T/Ar) ^ 1211n(T/AT) / ^ ' 

The only parameter A = kxjTf. is fitted according to the lattice data as A = 0.83 for model 
I and A = 0.7 for model II. The results of treatment (I) are shown in Figs. 2 and 3 for model 
I and model II, respectively. The inconsistency is transparently seen again. One can prove 
the deviation can not be conciliated by different choices of A. From Figs. 1-3, the failure of 
treatment (I) is obvious. The shapes of the curves are completely different. 

III. THE RESULTS OF TREATMENT (III) 



As demonstrated in our previous work [19"], the exact differential relations of thermody- 
namic functions in quasiparticle model should take the effective mass as a new independent 
degree of freedom. In grand canonical ensemble, the exact differential relation of f2 should 
be rewritten as Eqs.(3) and (4). We have proven that these formulae are self-consistent 
and agree with the statistical definitions and relations of thermodynamic quantities 
Employing Eq.(4), Eq.(8) and the energy definition 



g^mti = y ^ (14) 
I I e V » — i 

we calculate the corresponding thermodynamic quantities by treatment (III) for the SU(2) 
and SU(3) models in Sec. II. Our results are shown in Figs. 4-6. In Fig. 4, the solid curves and 
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the dashed curves refer to the lattice data and the results of treatment (III) for the SU(2) 
model, respectively. The solid curve and dashed curve for the entropy density coincide since 
the function m*(T) are determined from the lattice data of entropy 13j. Figs. 5 and 6 are the 
same as those of Figs. 2 and 3, except that the dashed curves are the results of our treatment 
(III). Comparing Figs. 1-3 from treatment (I) with the corresponding Figs. 4-6 from treatment 
(III), we see that the curves of our treatment (III) lies much closer to the lattice data, and 
the catastrophic deviation in Figs. 1-3 disappears. 

We are convinced by these results as well as Ref. 19j that treatment (I) with traditional 
partial derivative fails in the quasiparticle model with effective mass, and the thermody- 
namically consistent treatment (HI), by introducing the new independent variable m* in 
the partial derivative calculation, produces much better results. However, treatment (III) 
in its present form can not give some expected properties. We see from Figs. 5 and 6, the 
vanishing pressure at T = Tc is exhibited by the lattice data, since below the system 
becomes confined and no free particle exists any more. But the present result of pressure of 
treatment (III) for quasiparticle model is always finitely positive at T > 0. It is challenging 
to find a reasonable way to obtain a negative pressure modification to realize the phase 
transition at Tc within our consistent treatment, this will be the topic in next section. 



IV. MODIFICATION TO PRESSURE AND ENTROPY WITHOUT VACUUM 



Starting from the thermodynamic formula 



Q 

e + p - Ts = — = nfi = 0, 



(15) 



for qGP model and Eq.(4), we get 



e = T 



' dp ^ 
df 



— p = T 



(16) 



The difference between Eq.(16) and its usual form, e = — p, is significant. Instead 
of the total derivative we have the partial derivative with m* fixed, because m* is an 
independent variable in Eq.(3). We can integrate Eq.(16) and get 



p = T 



(17) 
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where f{m*) is an arbitrary integral function of m*. Eq.(17) shows the significant difference 
from the usual relation without taking m* as a free parameter. Obviously, if we choose 
f{m*) = 0, Eq.(17) reduces to the normal ideal gas case. 

According to Eq.(17), the thermodynamic potential fl, entropy S and X become 

rT 



Q ^ -pV^ -TV J ^dt - TVf{m*), 

+ f ^dt + f{m*) 



\ / V,m* 



e 
f 



(18) 
(19) 

(20) 



respectively. Then, noticing that — 0, the differential relation of internal energy U reads 

dU = d{n + TS) 

= TdS -pdV + Xdm* 



TV 



e 






f 






r 


1 


/ ds 






\ dm 



dV 



dt + TVf'{m*) 



Ty 



dm* 



e + Tj^^dt + Tf{m* 
+TV 



dV 



de ^ ,„ ^ ,m ^ f 



dm*dt + f\m*)dm* 



- Tf^dt + Tf{m*) 
sdV + Vds = d(Vs), 



dV -TV 



T,V 

T I ( de 



t"^ \ dm* 



dtdm* + f'{m*)dm* 



Ty 



(21) 



which means that U is not affected by the arbitrary function f{m*), and it still consists with 
the definition of internal energy in Eq.(14). This result is quite reasonable for the physical 
picture of the quasiparticle model. Prom Eqs. (17), (19) and (21), the relation 

o U + pV 



T 



(22) 



remains correct. 

Now we are in a position to determine f{m*) from the lattice simulation data. Hereafter 
we choose the SU(3) model I and model II as our working examples, taking advantage of 
the clearly expressed m* in Eqs. (11) and (12). 



With the effective mass in hand, the internal energy can be directly calculated from 
Eq.(14). For the pressure and entropy, one first calculate in equihbrium state as usual, the 
difference between the results and the lattice data gives the unknown function /(m*), as is 
shown in Eqs.(18) and (19). But since the internal energy obtained by the given m* may 
generally different from the lattice data, remembering U = TS — pV is satisfied by both 
the theoretical result and the data, the f{m*) fitted from the pressure and the one from the 
entropy will be different, as was shown in Figs. 7 and 8. Fig.7 refers to the SU(3) model I, 
where the dashed curve is for the f{m*)/T^ fitted from the entropy data and the dotted 
curve for the one fitted from the pressure data, respectively. Fig. 8 is the same as Fig.7 but 
for the SU(3) model II. We see from these two figures that the functions f{m*) obtained 
from the lattice pressure data and the entropy data lie quite close to each other. With these 
results, the modified pressure and entropy density for these two models are shown in Figs. 9 
and 10, with the dotted line representing the entropy density modified by f{m*) fitted from 
the pressure; and the dashed line representing the pressure modified by f{m*) fitted from 
the entropy, respectively. We see from these two figures that the curves of modified pressure 
and entropy density agree with the lattice data quite well. The internal energy density is 
not changed due to Eq.(21), so it remains the same as in Figs. 5 and 6 and is not plotted 
again in Figs. 9 and 10. It is worth pointing out that if m* is not given previously but fitted 
by the internal energy from Eq.(14), the f{m*) fitted from the pressure and the one from 
the entropy will be the same, and all the three curves will coincide exactly with the lattice 
data. In particular, the segments of negative pressure are plotted explicitly below Tc in 
Figs. 9 and 10. The function f{m*) contributes the negative pressure in a natural fashion. 
The appearance of the negative pressure, unlike treatment (II) jc, 13, 17], does not relies 
on the vacuum. This is an important character of our treatment (III). At last, we would 
like to emphasize again that the negative pressure shown in Figs. 9 and 10 at T < Tc does 
not means such states really exist. In fact the quasiparticle model is only applicable for 
T > Tc. We keep the segments of the curves below Tc just to demonstrate explicitly how 
the mechanism of negative pressure instability acts in the quasiparticle model within our 
treatment. 
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V. SUMMARY 



With the qGP model as a toy model, we demonstrate explicitly again the difficulty 
of treatment (I) with traditional partial derivative and prove that our treatment (III) is 
applicable, by the criterion of the lattice data. We also demonstrate the success of the 
newly introduced function f{m*) which can fit the lattice data self-consistently. It's worth 
emphasizing that the physical meaning of the additional function is completely different from 
the vacuum energy. This function does not appear in the expression of e, but change the 
values of p and s. The basic differences between our treatment (III) and other treatments 

n 

are as follows: (1). As was proven in Ref. [19|, treatment (III) is thermodynamically self- 
consistent and coincide with the statistics in equilibrium state. We do not need to introduce 
a vacuum and force it to cancel the thermodynamically inconsistent terms as in treatment 
(II). (2). We have an vacuum-unrelated integral function f{m*) in the consistent treatment, 
which can be fixed by fitting to experiment data or standard theoretical results. 



Together with our previous paper 191], we have established a self-consistent and complete 
thermodynamics for the effective mass quasiparticle model. Our treatment (III) is not only 
applicable to fermion system with nonzero chemical potential, but also to boson system with 
zero chemical potential. It can be wildly applied in all quasiparticle system. 
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FIG. 1: Internal energy density, pressure and entropy density as functions of temperature in SU(2) 
GP system. The solid curves are the lattice data from Ref.j25]; while the dotted curves are the 
results of treatment (I). 
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FIG. 2: Internal energy density, pressure and entropy density as functions of temperature in SU(3) 



GP system wit 
data from Ref. 



1 effective mass of model I and A = 0.83. The solid curves are the lattice simulation 



2g], and dashed curves are the results of treatment (I). 
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FIG. 3: The same as Fig.2, except that the effective mass is of model II, and A = 0.7. 
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FIG. 4: The same as Fig.l, except that the dashed curves are the results of treatment (III). 
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FIG. 5: The same as Fig.2, except that the dashed curves are the results of treatment (III). 
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FIG. 6: The same as Fig.3, except that the dashed curves are the results of treatment (III). 
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FIG. 7: The modification term f{m*)/T^ in SU(3) qGP model with m* of model I and A = 0.83. 
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FIG. 8: The same as Fig.7, except that the effective mass is of model II and A = 0.7. 
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FIG. 9: The same as Fig.5, but the pressure and entropy density are modified by f{m*) in Fig.7. 
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FIG. 10: The same as Fig.6, but the pressure and entropy density are modified by f{m*) in Fig.8. 
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